Resonance distribution in open quantum chaotic systems 
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In order to study the resonance spectra of chaotic cavities subject to some damping (which can 
be due to absorption or partial reflection at the boundaries), we use a model of damped quantum 
maps. In the high-frequency limit, the distribution of (quantum) decay rates is shown to cluster 
near a "typical" value, which is larger than the classical decay rate of the corresponding damped ray 
dynamics. The speed of this clustering may be quite slow, which could explain why it has not been 
detected in previous numerical data. 

PACS numbers: 05.45.Mt,03.65.Nk,42.55.Sa 



Recent experimental and theoretical studies have 
focussed on the dynamics of waves inside quasi 2- 
dimensional cavities which are "partially open"; this par- 
tial opening may be due to various physical phenomena. 
For instance, an acoustic wave evolving in air or in a 
metallic slab will lose intensity due to friction and heat- 
ing. In a microwave cavity, the dissipation mostly oc- 
curs at the boundary through ohmic losses. The light 
propagating inside a dielectric (micro) cavity is partially 
reflected at the boundary, which can be described as an 
"effective damping" at the boundary. In all these systems, 
the discrete stationary modes correspond to complex 
eigenvalues (or resonances) of the form k n = u> n — iT n /2, 
where r n is called the decay rate of the mode. 

When the shape of the cavity induces a chaotic ray dy- 
namics (e.g. the "stadium" shape), the eigenvalues {k n } 
cannot be computed analytically, but methods of "quan- 
tum chaos" can be applied to predict their statistical dis- 
tribution in the high-frequency limit uo n — > oo. Statis- 
tical studies of resonances started in the I960' with ini- 
tial applications to nuclear physics [l| . New applications 
emerged when experiments on mesoscopic quantum dots 
Q, microwave cavities 0| or optical fibers Q allowed to 
construct cavities with prescribed geometries, and study 
the dependence of the quantum dynamics with respect to 
this geometry. A recent interest for dielectric microcavi- 
ties comes from the potential applications to microlasers: 
choosing the shape of the cavity appropriately allows to 
produce a strongly directional emission [5j. The first step 
to understand the (nonlinear) lasing modes is to study 
the passive (resonant) modes of the cavity. 

Various dissipation effects have been taken into ac- 
count by adding to the self-adjoint Hamiltonian (repre- 
senting the dissipationless system) an effective imaginary 
part, which describes the coupling between the internal 
cavity modes and the external channels [(| . One analyti- 
cal tool to study chaotic cavities has been to replace the 
Hamiltonian (and sometimes also the effective coupling) 
by some sort of random matrix: this has lead to theoret- 
ical distributions, which have been favorably compared 
with numerical or experimental spectra [7|, |8j . 

In this Rapid Communication we focus on situations 
where the coupling is strongly nonperturbative, and is 



distributed over a large part of the cavity or of its bound- 
ary, so that the number of coupled channels becomes 
"macroscopic" in the high-frequency (semiclassical) limit. 
Using a nonrandom model of damped quantum maps, we 
find that, in this limit, the distribution of quantum decay 
rates becomes asymptotically peaked on a "typical" value 
Jtyp, which is the ergodic mean of the local damping rate. 
This clustering does not seem to appear if one replaces 
the unitary part of the quantum map by a random uni- 
tary matrix, as is often done in the quantum chaos lit- 
erature [13, [2a|- Such a clustering has been ri goro usly 
proved for damped waves on ergodic manifolds [1J|; we 
believe it to occur as well in the various types of partially 
open quantum systems mentioned above. Yet, the width 
of the distribution may decay very slowly in the semi- 
classical limit (< (ln/c) -1 / 2 ), which could explain why 
this semiclassical clustering is hardly visible in numerical 
computations of chaotic dielectric cavities 0, [3, E2] or 
damped quantum maps fla |: such a slow decay indeed 
occurs within a solvable toy model we briefly describe at 
the end of this note. 

Let us now describe the model of damped quantum 
map, which has been introduced and numerically in- 
vestigated in [HI to mimick the resonance spectra of 
dieletric microcavities. To motivate this model, we first 
briefly analyse the dynamics of a few cavity wave sys- 
tems. The first situation consists in a smooth absorption 
inside the cavity, represented by the damped wave equa- 
tion (<9 t 2 - A + 2b{x)d t )ip{x,t) = 0. Here the "damp- 
ing function" b{x) > measures the local absorption 
rate. A high-frequency wavepacket evolving along a 
classical trajectory is continuously damped by a factor 
~ e - So b( x ( s )) ds . The classical limit of the dynamics con- 
sists in the propagation of rays with decreasing intensity, 
also called weighted ray dynamics (Fig. [H (a)). When 
the dissipation occurs at the boundary (e.g. through 
ohmic losses) , an incident high-frequency wavepacket hit- 
ting the boundary will be reflected, with its amplitude 
reduced by a subunitary factor a(q,(p) (Fig.[lJ (b)). The 
same phenomenon effectively occurs in the case of light 
scattering through a quasi-2D dielectric microcavity of 
optical index n > 1. The rays propagating inside the 
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FIG. 1: (Color online) (a) Weighted ray dynamics inside a 
cavity with inhomogeneous absorption, (b) absorption (or 
partial reflection) at the boundary. The ray intensity corre- 
sponds to its thickness. Dashed lines correspond to refracted 
rays. 



cavity are partially reflected at the boundary, the re- 
maining part being refracted outside and never return 
provided the cavity is convex (dashed lines in Fig. [TJ 
(b)). In the high-frequency limit, the reflection factor is 
given by Fresnel's coefficient, which depends on the light 
polarization and of p = sin <p. For instance, in case of 
transverse magnetic polarization the coefficient is a sim- 
ple complex function citm(p) [13j, which has unit mod- 
ulus when \p\ > 1/n (full reflection) and is minimal at 
«tm(0) = 

To analyze a 2D classical billiard, it is convenient to 
reduce the flow to the bounce map n : (q,p = sin^) i— » 
(q',p' — sm(p'), which acts canonically on the bound- 
ary phase space. At the quantum level, the spectrum 
of the closed cavity can be obtained by studying a k- 
dependent integral operator acting on the boundary, 
which effectively quantizes the bounce map, with an ef- 
fective Planck's constant K e ff = k^ 1 [1 51 ] - 

This observation leads to consider canonical maps k 
on simple 2-dimensional phase spaces, and quantize them 
into unitary propagators (quantum maps) Un{k) of finite 
dimension N ~ h~ff A Gaussian wavepacket \q,p) 

localized at the phase space point (q,p) is first trans- 
formed unitarily into a deformed wavepacket Un\q,p), 
localized near n(q,p). To induce some damping, we 
then multiply this state by a factor a(n(q,p)), which can 
be implemented by applying to Un \q,p) the operator a 
quantizing the damping factor. The latter is generally 
complex- valued, we will assume that it satisfies the fol- 
lowing bounds: 



Vq,p, < a min < \a(q,p)\ < a max = 1 



(1) 



These two steps lead to the definition of the damped quan- 
tum map 

Mn = Mjy(a, k) — do Un(k) . 

The classical limit of the dynamics generated by Mjy acts 
on "weighted point particles": a point at position (q,p) 
is moved to n(q,p) and its weight is reduced by a fac- 
tor \a(K(q,p))\ 2 . This is the discrete-time version of a 
weighted ray dynamics. Compared with cavity systems, 



this model has two main advantages: one can easily en- 
gineer a map k with specific dynamical properties; the 
spectrum of M N (a, n) is easier to study both numerically 
and analytically. 

The spectrum {^-j}i<j<N of Mn(ci, k) is the main 
object of our study (eigenvalues are ordered by decreasing 
moduli). To compare it with the resonance spectrum 
of a damped cavity, one should extract from the latter 
an interval {\u n — k\ < tt} around the frequency k ~ 
N. The distribution of the decay rates {r„ : \u/ n — 
k\ < 7r} should parallel that of the decay rates = 

-21n|Af ] \}i<j<n. 

A similar model had been introduced in [3, EH] to 
mimick "fully open" cavities: the damping factor a(z) 
was then vanishing inside the opening. Such systems 



were characterized by a fractal Weyl law [121, LL8|, Ll9fl: the 
number of resonances in a strip {\uj n \ < k, T n < T} grew 
like k 1+s , where 6 < 1 was given by the fractal dimen- 
sion of the trapped set. On the opposite, the bounds 
([I]) imply that for N large enough, Afjv is invertible, 
and its N eigenvalues are contained inside the annulus 
{a m in < \^j N ^\ < !}• Transposed to the case of an 
absorbing cavity, it implies that all high-frequency res- 
onances are contained in a fixed strip {r„ < r max }, and 
the number of modes N{k n : \u) n \ < k} asymptot- 
ically grows like C k 2 , thus satisfying a standard Weyl 
law [IJ]. The situation is more complicated for dielec- 
tric cavities. Explicit solutions in the case of the circu- 
lar cavity [ll| suggest that resonances split between two 
well-separated groups: "inner" resonances contained in a 

1 /3 

strip {r„ < T}, and "outer" resonances T n ~ luJ as- 
sociated with modes localized outside the cavity. Since 
our damped quantum map only acts on states localized 
"inside the torus", we believe that the above Weyl asymp- 
totics correctly counts the inner resonances of dielectric 
cavities (the fractal Weyl law recently observed in [l2j is 
probably a finite-frequency artifact). 

To obtain a more precise description, one needs to it- 
erate the dynamics, that is study the time-n evolution 
Mj!f: Applying the quantum-classical correspondence 
("Egorov's theorem"), we find that 



(M N (a, K pM N (a, K ) n ) 1/2n n a 



(2) 



where the function a n = (IliLi l a n% 

\) 1/n is the aver- 
age damping over trajectory stretches of length n. The 
approximation is valid in the semiclassical limit N — > 00. 

Much can be drawn from the knowledge of the func- 
tions — 21na„ in the long-time limit n 3> 1. Their ranges 
consist in intervals I n (a) C 7„_i(a), which converge to a 
limit interval loo (a) when n — > 00. The above identity 
implies that the quantum decay rates 7^ must be con- 
tained in /00(a) for large enough N [Tiff . Numerical [l3| 
and analytical [20] studies indicate that the "quantum 
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ranges" Jjv(a) = [yj , J^] generally remain strictly 
inside /00(a), in particular stay at finite distance from 
zero. Adapting methods used to study scattering sys- 
tems 22, 123J , one finds that high-frequency decay rates 

should be larger than j gap == f — 2"P K (m \a\ — A"/2), where 
V K {-) is the topological pressure associated with the map 
k and the observable (In \a\ - X u /2) [H], \ u (z) being the 
expansion rate of K along the unstable direction. The 
lower bound j gap may be trivial (negative) when \a(z)\ 
varies little across the phase space (see the last line in 
Table U). 

Since K is chaotic, the value distribution of 
— 21na„(z) becomes peaked around its average jt yp — 
—2 Jin \a(z)\ dz when n — » 00. The central limit theo- 
rem [25| shows that this distribution is asymptotically a 
Gaussian of width around Jtyp- This distribution is 
semiclassically connected with the spectral density of a n : 
denoting {s\ N }i<i<N its eigenvalues and V the volume 
on the torus, we have the Weyl law 



< S }«iVV(a- 1 ([o 



»])) 



(3) 



From the LHS approximately counts the singular 
values of the matrix Mjy. Using the Weyl inequalities 

[iij ]. we obtain that most of the decay rates {jj }i<j<N 

satisfy 7^ > j typ - (■ 

Applying the same argument to the inverse quantum 
map M^ 1 M7v(a~ 1 °K, k -1 ), we eventually find that in 
the semiclassical limit, most decay rates cluster around 
Jtyp, which we thus call the typical decay rate. More 
precisely, the fraction of the decay rates {7^ } which 
are not in the interval [jtyp — £, Jtyp + s] goes to zero 
when N — * 00. 

By pushing the quantum-classical correspondence up 
to its limit, namely the Ehrenfest time n ~ Cm AT, we 
find that the width of the decay rate distribution is at 
most of order (In AQ -1 / 2 (a rigorous proof will be given in 
[13]). Our numerics (see Fig[3]) are compatible with this 
upper bound. Such a slow decay could explain why this 
concentration has not been detected in previous studies. 
For a solvable toy model presented at the end of this note, 
the distribution will be shown to be indeed a Gaussian 
of width - CQnN)- 1 / 2 . 

Let us compare the quantum decay rates with the "clas- 
sical decay rate" j c i of the corresponding weighted dy- 
namics. The latter, introduced in [ilj in the framework 
of dielectric microcavities, is obtained by evolving an ini- 
tial smooth distribution of points through the weighted 
dynamics: for large times n, the total weight of the distri- 
bution decays like W e - " 7ci . As in the case of fully open 
systems [13], j c \ can be expressed as the topological pres- 
sure j c i = —V K (2 In \a\ — A"). Convexity properties of the 
pressure allow to compare this classical decay rate with 
the two rates obtained above: j gap < Jd < Jtyp, and the 
inequalities are generally strict (see table HJ. The quan- 
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FIG. 2: (Color online) Spectra of Mjv(k,o<) for N = 2100 
and damping factors a\ (a) and 02 (b). We plot {iln(\^ N ^)} 
to mimick the spectrum of a damped cavity near k ^ N (the 



vertical coordinates correspond to 
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/2). The horizontal 



lines indicate — 7*/2 for the theoretical rates given in tabled] 
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7ci 


Jtyp 


Oi 


0.715 


1.079 


2.774 


ffll 


0.734 


1.079 


3.070 


ffl2 


-0.523 


0.521 


0.633 



TABLE I: Values of the theoretical rates for k and the various 
damping functions we use. 



turn ranges Jn{c) may or may not contain the classical 
rate j c i (see Fig. [2]). 

The map we consider in our numerics is the 3-baker's 
map, which acts canonically on the 2-dimensional torus 
{(<l,p) £ [0, l) 2 }- It is given by K,{q,p) — (3g mod 
1, p+ | 3g ] ), and generates a strongly chaotic dynamics. 
This map is quantized as in [28] , into a sequence of uni- 

Gn/3 



tary matrices Un = G 



N 



G 



N/3 



where 



G 



N/3/ 



(G M ) jk = exp ( - + l)(k + I)) is the "sym- 

metrized" discrete Fourier transform. We choose damp- 
ing factors of the form a(q), so their quantizations a are 
simply diagonal matrices with entries a((j + l/2)/N). 
The factor a\{q) has a plateau a\{q) = 1 for q G 
[1/3,2/3], another one ai(q) = 0.1 for q £ [0,1/6] U 
[5/6,1], and varies smoothly inbetween. It approxi- 
mates the piecewise constant function a\{q) which takes 
values 0.1, 1, and 0.1 respectively on the intervals 
[0,1/3), [1/3,2/3) and [2/3,1). Our second choice is the 
"smoother" function a 2 (g) = 1 — sin(27rg) 2 /2. Since we 
use a single map k, the damped quantum maps will be 
abbreviated by Mpj(ai). 

We first notice that all these factors reach their ex- 
tremal values a m i n , a max on the fixed points (0,0) and 
(1/2, 1/2) of k. As a result, for each of them the asymp- 
totic range /00(a) is equal to [a m i n , a max ]. The theoretical 
rates j ga p, Jci and jt yp for these three factors are given 
in Table HI In Fig [2] we plot the spectra of MAr(aj) for 
N = 2100 (the theoretical bound j gap for 02 is negative, 
hence irrelevant). We check that all quantum rates are 
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larger than 7 9a p- In the case of Mry(ai), all quantum 
rates are also larger than j c i, while M^^) admits a few 
smaller decay rates. 

The clustering of decay rates around j tyP is already 
perceptible in FigO To make it more quantitative, in 
FiglU (a) we plot the cumulative distributions of decay 
rates. At first glance, the widths of the distributions 
around jt VP seem to depend little on N. Enlarging the 
set of data, we plot these widths on Fig[3] (b). They 
indeed decay with AT. The 2-parameter power-law fits 
A' N~ B lead to small exponents B' , which seem to favor 
the logarithmic fits A (lnN)~ B : the latter decay slightly 
faster than the theoretical upper bound (In TV) -1 / 2 . 
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FIG. 3: (Color online) (a) Cumulative decay rate distribution 
for Mjv(k, ai). The vertical bars indicate the rates ytyp- (b) 
Widths of the decay rate distributions, together with the best 
2-parameter fits A(lniV) -s and the asymptotic standard er- 
rors. The power-law fits A' N~ B ' yield A' = 1.35 ± 5% , B' = 
0.09±8% for Mjv(oi) and A' = 0.2±18%,B' = 0.10±26.5% 
for M N (a 2 ). 



It is possible to construct a solvable quantization of the 
baker's map by taking the quantum parameter N = 3 fc , 
k € N, and replacing the discrete Fourier transform 
Gjv by the Walsh transform If we then select 

a damping factor which, like a\, takes constant values 
a J on the intervals q G [j/3, (j + l)/3), the quantum 
model remains solvable. The spectrum of Mm relies on 
the eigenvalues {Ai} of the 3x3 matrix diag(a J )G ! J . 
Taking 7, = — 21n|Aj|, the N quantum decay rates 
can be indexed by the sequences 77 = nxTh. . . . rjk, with 
t]i £ {1,2,3}: they are given by ^ = £ Em=i 7»? m ■ For 
instance, in the case of the damping function 5i, the 
rates ji take the values (0.803, 3.801, 4.605). From this 
explicit expression, one easily draws that, in the limit 
k — > 00, the distribution of the {7^} converges to a 
Gaussian of average jt yp — (IZi=i7i)/^ an< i variance 
AE-=i(7*-7^) 2 = C(mA0-i. 

To summarize, we have studied the spectra and eigen- 
modes of damped quantum chaotic maps, a toy model for 
various types of partially open quantized chaotic cavities, 
in a regime where the damping is both macroscopic and 
strongly nonperturbative. We have shown that the quan- 



tum decay rates remain inside a fixed interval, and that 
most of them cluster around the mean damping rate Jtyp- 
These statistical properties seem to differ from those of 
non-hermitian random matrices used to represent such 
open systems. 
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